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■ Abstract: "Fuzzy CP^", which is a four-dimensional fuzzy manifold analogous to the 
fuzzy 2-sphere (S^), appears as a classical solution in the dimensionally reduced 8d Yang- 

^ . Mills model with a cubic term involving the structure constant of the SU(3) Lie algebra. 

Although the fuzzy S^, which is also a classical solution of the same model, has actually 
smaller free energy than the fuzzy CP^, Monte Carlo simulation shows that the fuzzy 
CP2 is stable even nonperturbatively due to the suppression of tunneling effects at large 
N as far as the coefficient of the cubic term (a) is sufficiently large. As a is decreased, 
both the fuzzy CP^ and the fuzzy collapse to a solid ball and the system is essentially 
described by the pure Yang-Mills model (a = 0). The corresponding transitions are of first 
order. The gauge group generated dynamically above the critical point turns out to be of 
rank one for both CP^ and cases. Above the critical point, we also perform perturbative 
calculations for various quantities to all orders, taking advantage of the one-loop saturation 
of the effective action in the large-A^ limit. By extrapolating our Monte Carlo results to 
N = oo, we find excellent agreement with the all order results. 
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1. Introduction 

Fuzzy spheres [1], which are simple compact noncommutative manifolds, have been recently 
discussed extensively in the literature. One of the motivations comes from the general 
expectation that noncommutative geometry provides a crucial link to string theory and 
quantum gravity. Indeed Yang-Mills theories on noncommutative geometry appear in a 
certain low energy limit of string theory [2]. There is also an independent observation 
that the space-time uncertainty relation, which is naturally realized by noncommutative 
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geometry, can be derived from some general assumptions on the underlying theory of 
quantum gravity [3] . Another motivation is to use fuzzy spheres as a regularization scheme 
alternative to the lattice regularization [4]. Unlike the lattice, fuzzy spheres preserve the 
continuous symmetries of the space-time considered, and hence it is expected that the 
situation concerning chiral symmetry [5-23] and supersymmetry might be improved. 

As expected from the connection to string theory [24], fuzzy spheres appear as classical 
solutions in matrix models with a Chern-Simons-like term [25-29] and their dynamical 
properties have been studied in refs. [30-39]. One can actually use matrix models to 
define a regularized field theory on the fuzzy spheres as well as on a noncommutative 
torus [40] , which enables nonperturbative studies of such theories from first principles [41] . 
These matrix models belong to the class of the so-called dimensionally reduced models (or 
large- iV reduced models), which is widely believed to provide a constructive definition of 
superstring and M theories [42-44]. The space-time is represented by the eigenvalues of the 
bosonic matrices, and in the IIB matrix model [43], in particular, the dynamical generation 
of /our-dimensional space-time (in ien-dimensional type IIB superstring theory) has been 
discussed by many authors [45-58]. 

In rcf. [59] we have studied the dimensionally reduced 3d Yang-Mills-Chern-Simons 
(YMCS) model, which has the fuzzy 2-sphere (S^) as a classical solution [26]. Unlike 
previous works we have performed nonperturbative first-principle studies by Monte Carlo 
simulation. We observed a first-order phase transition as we vary the coefficient (a) of the 
Chern-Simons term. For small a the large- A?^ behavior of the model is the same as in the 
pure Yang-Mills model, whereas for large a a single fuzzy appears dynamically. 

For obvious reasons it is interesting to extend this work to a matrix model which 
accommodates a /owr-dimensional fuzzy manifold. In ref. [60] we have studied the dimen- 
sionally reduced 5d Yang- Mills model with the quintic Chern-Simons term [27, 29] , which 
is known to have the fuzzy 4-sphere (S^) as a classical solution [61]. Unlike the fuzzy 
case, however, the fuzzy is unstable at the classical level, and Monte Carlo simulation 
confirmed that it does not stabilize even at the quantum level. The negative result is es- 
sentially due to the fact that the quintic Chern-Simons term has higher powers in Afj, than 
the Yang-Mills term. 

This motivates us to return to the class of models with a cubic term. As a candidate of 
a four-dimensional fuzzy manifold, we study the fuzzy CP^ [28,34,62-65], which appears 
as a classical solution in the dimensionally reduced 8d Yang-Mills model with a cubic term 
involving the structure constant of the SU(3) Lie algebra. In fact the fuzzy S^, which is also 
a classical solution of this model, has smaller free energy than the fuzzy CP^. Monte Carlo 
simulation shows, however, that the fuzzy CP^ is stable even nonperturbatively due to the 
suppression of tunneling effects at large N as far as the coefficient (a) of the cubic term 
is sufficiently large. As we decrease a, both the fuzzy CP^ and the fuzzy S^ collapse to a 
solid ball, and the system is essentially described by the pure Yang- Mills model (a = 0). 
The corresponding phase transitions are of first order, and the lower critical point agrees 
with the analytical result obtained from the one-loop effective action. Since the one-loop 
effective action is saturated at one loop in the large- limit, the analytical result for the 
critical point is expected to be free from higher loop corrections. 
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Above the critical point, we perform perturbation calculations for various quantities to 
all orders, taking advantage of the one-loop saturation of the effective action in the large- iV 
limit. This technique was originally proposed for a supersymmetric model [37], where the 
effective action is saturated only at two loop. In the bosonic case, one can obtain all order 
results by performing essentially the one-loop calculation [66] . By extrapolating our Monte 
Carlo results for various observables to iV = oo, we find excellent agreement with the all 
order results. 

In the large- iV reduced models, not only the space-time [45] but also the gauge group 
[67] is expected to appear dynamically. While there are certain evidences in the IIB matrix 
model that indeed four- dimensional space-time appears dynamically [53-55,57,68,69], the 
issue of the gauge group is totally unclear. The models we are studying may be considered 
as a toy model in which one may obtain a definite answer to such a question, since the 
gauge group of rank k naturally appears if the true vacuum is given by k coincident fuzzy 
manifolds. The value of k should be determined dynamically, and the result may of course 
depend on the model one considers. In the present model we find the gauge group to be 
of rank one for both the fuzzy CP^ and the fuzzy S^. In arriving at this conclusion, the 
existence of the first-order phase transition plays a crucial role as in our previoTis work [59] . 

This paper is organized as follows. In section 2 we define the model and show that the 
fuzzy CP^ and the fuzzy appear as classical solutions. In sections 3 and 4 we study the 
properties of the fuzzy CP^ and the fuzzy S^, respectively. In section 5 we compare the free 
energy for the fuzzy CP^ and the fuzzy to discuss which is the true vacuum whenever 
they exist. In section 6 we study the properties of the k coincident fuzzy CP^. In section 
7 we determine the rank of the dynamical gauge group for the fuzzy CP^ and the fuzzy 
S^. Section 8 is devoted to a summary and discussions. In appendix A we present the 
explicit form of the fuzzy CP^ configuration. In appendix B we formulate the perturbation 
theory around fuzzy manifolds and obtain an expression for the one-loop free energy. In 
appendices C and D we show the perturbative calculations for the fuzzy CP^ and the fuzzy 
S^, respectively. 



2. The model and its classical solutions 

The model we study is defined by the action 

S = Nti:(^~[A^,A^]'^ + ^iaf^^pA^A^Ap^ , (2.1) 

where A^ {n = 1, - ■ ■ ,8) are N x N traceless hermitian matrices. Here and henceforth we 
sum over repeated indices. The coefficient /^^p is the structure constant of the SU(3) Lie 
algebra, whose nonzero components are given explicitly by 

\/3 1 

/l23 = 1 , /458 = /678 = ' /l47 = /246 = /257 = fub = /5I6 = fm7 = ■ (2-2) 

The pure Yang-Mills model (a = 0) and its obvious generalization to D dimensions 
with D matrices A^ = 1, . . . , D) have been studied by many authors. In particular, the 
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large- dynamics of the model have been studied by the 1 / D expansion and Monte Carlo 
simulation [70]. The partition function was conjectured [71] and proved [72] to be finite for 
N > D/ {D — 2). (See refs. [73-75] for the super symmetric case.) The partition function in 
the presence of the Chern-Simons term has been studied analytically for iV = 2 [76], and 
it turned out to be convergent in the supersymmetric case, but not in the bosonic case. 
It is also proved that adding a Myers term (the cubic term in the present case) does not 
affect the convergence as far as the original path integral converges absolutely [77], which 
means, in particular, that the partition function of our model is convergent for > 4. 
The classical equation of motion of the model (2.1) is given by 

[A,, [A^, A,]] - i a U,p[A„ Ap] = 0. (2.3) 

One can easily see that there exists a solution of the form 

Ai, = aTi,, (2.4) 

where satisfies the SU(3) Lie algebra 

[T^,n]=if^,pTp. (2.5) 

Hence one obtains a classical solution for each of the AT-dimensional representations of the 
SU(3) Lie algebra. Let us consider the case in which is given by the irreducible (m, n) 
representation Tjf^'"'^ . Such a solution exists when the size of the matrices A^ is 

iV = ^(m+l)(n + l)(m + n + 2) . (2.6) 

The explicit form of J"^™'") jg given in the appendix A. 

The space represented by the matrices Af^ = aT^™'"'' has SU(3) isometry. There 
are two kinds of manifold whose isometry is SU(3). One is SU(3)/U(2), and the other 
is SU(3)/(U(1) X U(l)). In fact the SU(3)/U(2) space corresponds to the (m,0) or the 
(0, n) representation, whereas the SU(3)/(U(1) x U(l)) space corresponds to the (m, m) 
representation. 

In what follows we consider the CP^ = SU(3)/U(2) space, which is given by 

4CP') = ar^-'°) . (2.7) 

Note that this solution exists only for 

iV = i(m + l)(m + 2) = 3,6, 10, 15,21,--- . (2.8) 

This is in contrast with the fuzzy S^ case [59] , where the corresponding solution exists for 
arbitrary N. As in ref. [59] we define the "radius-squared matrix" 

Q = {A^? , (2.9) 
which is useful for distinguishing various solutions. Using the identity (A. 8), we obtain 



Q = p'lN, P = a\M^ (2.10) 
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for the CP^ solution, which implies that the fuzzy CP^ has the radius p. 

We note that the model also has the fuzzy type of classical solutions. As an example, 
let us consider 



^(S2) = J " ^M^^ for = 1, 2, 3 , 
I otherwise , 



(2.11) 



where L^^^ is the iV-dimensional irreducible representation of the SU(2) Lie algebra. The 
radius-squared matrix Q is given in this case as 



(2.12) 



Although there are other fuzzy solutions with a smaller radius, the one given above is 
considered the most relevant since it has the smallest action among them. 

3. Properties of the fuzzy CP^ 

In order to study the properties of the fuzzy CP^, we simulate the model (2.1) using (2.7) 
as the initial configuration. Wc apply the heat bath algorithm developed in ref. [70] for 
the pure Yang-Mills model by implementing the cubic term as explained in the appendix 
A of ref. [59]. 
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Figure 1: The quantities -^{S) and -^(■i-tr(A^)^) are plotted against a = aN^ for N = 
10,15,21,28 (m = 3,4,5,6) for the fuzzy CP^ start. The dotted, dashed and solid lines represent 
the classical, one-loop and all order results, respectively, at large A''. 

From perturbative calculations, it turns out natural to fix the rescaled parameter 

a = am , (3.1) 

when wc take the large- limit. In figure 1 we plot the Monte Carlo results for the action 
{S) and the "extent of space-time" (■^tr(74^)^) (with an appropriate normalization factor) 
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Figure 2: The quantities -^{S) (left) and -^{j^tr{Ai^)'^) (right) are plotted against for a = 3.0 
for the fuzzy CP^ start. The dotted lines represent a linear fit, and the horizontal solid (dashed) 
lines represent the all order (one- loop) results obtained analytically at iV = oo. The large- iV 
extrapolation demonstrates excellent agreement with the all order result. 

against the rescaled parameter a ioi N = 10,15,21,28 (m = 3,4,5,6). We observe a 
discontinuity around 

a = 4?^'^ ^ 2.3 , (3.2) 

which suggests the existence of a first-order phase transition. An analogous first-order 
phase transition has been found also in the 3d YMCS model [59]. The critical point (3.2) 
agrees well with the analytical result (C.IO) with k = I. 

Below the critical point, the Monte Carlo results are almost independent of a. This is 
because the cubic term in the action (2.1) takes small values, and hence it does not play any 
role in this regime. Note, in particular, that ;^(tr(A^)^) ~ 0(1) as in the pure Yang-Mills 
model (a = 0) [70], which means that j^{tT{Aij,)'^) has different large- AT behaviors in the 
two phases, as clearly seen in figure 1. 

Above the critical point, we observe that Monte Carlo results approach the all order 
results (C.40), (C.41) with fc = 1 as increases. In order to clarify the finite- effects, in 
figure 2 we plot the same quantities against for = 10, 15,21,28 (m = 3,4,5,6) with 

fixed a = 3.0 > acr • The data can be nicely fitted to a straight line, which implies that 
the leading finite- AT effect is of 0(;^). This allows us to make a reliable extrapolation to 
N = oo, and we find excellent agreement with the all order results. 

4. Properties of the fuzzy 

In this section we study the properties of the fuzzy by Monte Carlo simulation us- 
ing (2.11) as the initial configuration. Perturbative calculations suggest that the natural 
definition of the rescaled parameter in this case is 

a = aN^ , (4.1) 
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Figure 3: The quantities -^{S) and -i- (-i-tr(Aju)^) are plotted against a = aN^ for N = 
10,15,21,28 with the fuzzy start. The dotted, dashed and soHd Hues represent the classical, 
one-loop and all order results, respectively, at large A'". 
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Figure 4: The quantities ^(5) (left) and -i-(^tr(A^)^) (right) arc plotted against for a = 5.0 
with the fuzzy start. The dotted lines represent a linear fit, and the horizontal solid (dashed) 
lines represent the all order (one- loop) results obtained analytically at A'' = oo. The large- A'' 
extrapolation demonstrates excellent agreement with the all order result. 



unlike (3.1) in the fuzzy CP^ case. In figure 3 we plot -^{S) and ;^(;^tr(A^)^) against a. 
We observe a discontinuity at 

a = a^f ) ~ 3.2 , (4.2) 



which suggests the existence of a first-order phase transition. The critical point (4.2) agrees 

well with the analytical result (D.6) with k = 1. In terms of the unrescaled parameter a, 
the critical point for the fuzzy CP^ and the fuzzy are Oc?^ ^ — ^^ifi and ■* ~ -^^j 

respectively, which means that Oci < Ocr • Below the critical point, Monte Carlo 
results are identical to those for the fuzzy CP^ start presented in the previous section, as 
expected. 
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Above the critical point a > cxcy , Monte Carlo results arc quite close to the all order 

results (D.13), (D.14) with k = 1. In figure 4 we plot the two quantities against for 

(S'^) 

N = 6, 10, 15,21 with fixed d = 5.0 > dcr ■ Monte Carlo results can be nicely fitted to 
a straight line, which implies that the leading finite- A?" effect is 0(-^), as compared with 
0(;^) for the fuzzy CP^ case. The large- AT extrapolation demonstrates perfect agreement 
with the all order results obtained in the large- A" limit. 

5. CP^ versus — which is the true vacuum? — 

In the previous two sections, we have seen that both the fuzzy CP^ and the fuzzy are 
stable for sufficiently large a. In this section we discuss which of the two describes the true 
vacuum. For that purpose we compare the free energy for the fuzzy CP^ and the fuzzy S^, 
which are obtained to all orders in perturbation theory in the large-A limit as (C.37) and 
(D.12), respectively. Setting k = 1 and rewriting in terms of the unrescaled parameter a, 
the free energy reads 

^(cp^) ^ n2(_^^q,^^ lo^m') -— - ^ (5 1) 

~ \ 6 ^Dioga + iog^oiv ; ^^^^ 2a^^m J'^ ' 

' " ^ (-^ + 61oga + logAr- - -^--^--^-...j. (5.2) 

Thus in the region a > ac?^ ^ = ^^iyi , where both the fuzzy CP^ and the fuzzy exist, 
we find that the fuzzy has smaller free energy than the fuzzy CP^ at large A". Therefore 
the fuzzy CP^ cannot be the true vacuum in the present model. 

We note, however, that the fuzzy CP^ appears to be very stable. For instance, we have 
performed a simulation with the fuzzy CP^ start for A = 10 (m, = 3) and a = 1.4, which 
is just above the critical point. The fuzzy CP^ does not decay into the fuzzy even after 
10^ sweeps of the heat bath algorithm. This suggests the existence of a potential barrier 
between the two vacua, which presumably increases with N. We therefore consider that 
the fuzzy CP^ stabilizes due to the suppression of tunneling effects in the large- AT limit. 

6. Properties of the k coincident fuzzy CP^ 

In this section we discuss the properties of the k coincident fuzzy CP^ configuration 

4'=CP^)^ar(--°)0l,, (6.1) 
which is also a classical solution of the model. The size of the matrices should now be 

N = ^k{m + l){m + 2) . (6.2) 

Such a configuration is important since it gives rise to a gauge theory on the fuzzy CP^ 
with the gauge group of rank k. 

We have performed Monte Carlo simulation with the initial configuration given by 
(6.1) with k = 2. In figure 5 we plot the action S and the eigenvalues of the radius-squared 
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Figure 5: The history of the action S (left) and the eigenvalues of the radius-squared matrix Q 
(right) for a = 2.0 and = 42 using the k = 2 coincident fuzzy CP^ as the initial configuration. 
The horizontal line in the left (right) plot represents the one-loop (classical) result for the k = 2 
coincident fuzzy CP^ at a = 2.0 and N = 42. 
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Figure 6: The quantities -^{S) and "^^("^tr(A^)^) are plotted against a = aN^ for N = 
20,30,42 (m = 3,4,5) with the k = 2 coincident fuzzy CP^ start. The dotted, dashed and sohd 
lines represent the classical, one-loop and all order results, respectively, at large N. 



matrix Q defined by (2.9) against the number of "sweeps" in the heat bath algorithm [70]. 
We observe that the k = 2 fuzzy CP^ decays after 600 sweeps. 

Although the k = 2 coincident fuzzy CP^ is thus only meta-stable, we may measure 
various observables before it actually decays in Monte Carlo simulation. In figure 6 we plot 
the results against a for iV = 20, 30, 42 (m = 3, 4, 5). We observe a discontinuity at 

a = «(^=2CP^)^2.7, (6.3) 

which agrees with the analytical result (C.IO). Above the critical point, our Monte Carlo 
results agree well with the all order results (C.40), (C.41) with k = 2. 



-9- 



7. Dynamical gauge group 

In the previous section, we have seen in Monte Carlo simulation that the k coincident fuzzy 
CP^ is unstable in the k = 2 case. This instability is related to the zero modes that appear 
in the perturbation theory around these configurations. (See appendix C.l.) Although we 
cannot exclude the possibility that the instability disappears in the large-iV limit, we show 
that the multi-fuzzy CP^ cannot be the true vacuum anyway, by comparing the free energy 
calculated omitting the zero modes. The explicit form of the free energy to all orders in 
perturbation theory is obtained as (C.37) above the critical point (C.IO), which we plot for 
A; = 1, • • • , 6 in figure 7 (left). We find that the k = 1 case gives the smallest free energy 
for all values of a. Thus, we conclude that the dynamical gauge group for the fuzzy CP^ 
is of rank one. 

We repeat the same analysis for the fuzzy case. The free energy for the k coincident 

fuzzy is given to all orders in perturbation theory above the critical point (D.6) by 
(D.12), which wc plot for k = 1, • • • ,6 in figure 7 (right). We find that the k = 1 case 
gives the smallest free energy for all values of a. Therefore, the dynamical gauge group is 
of rank one in this case as well. 

Let us emphasize that the existence of the first-order phase transition plays an im- 
portant role in arriving at these conclusions. As one can see from the free energy (C.37), 
(D.12), the classical term favors small k, while the one-loop term proportional to log A; fa- 
vors large k. Therefore, if we disregard the existence of the first-order phase transition, the 
rank of the dynamical gauge group could be A; > 1 at small a. What actually happens is 
that the critical point (C.IO), (D.6) increases with k, and since we have to restrict ourselves 
to the region above the critical point, the one-loop term cannot compensate the effect of 
the classical term. 




Figure 7 : The free energy obtained to all orders in perturbation theory is plotted above the critical 
point for the k coincident fuzzy CP^ (left) and fuzzy (right) with fc = 1, • • • ,6. We have taken 
the large- iV limit after subtracting the irrelevant constant term proportional to \ogN. In both 
cases, the k = 1 case represented by the solid line gives the smallest free energy. 
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8. Summciry and discussions 

In this paper we have apphed nonperturbative techniques and ideas, which have been 
developed in ref. [59], to a four-dimcnsional fuzzy manifold, the fuzzy CP^. The present 
model may be considered as a natural extension of our previous model in the sense that 
the epsilon tensor, which is nothing but the structure constant of SU(2), is now replaced 
by that of SU(3). Since the number of bosonic matrices should be equal to or larger 
than the dimensionality of the algebra, it is taken to be 8 instead of 3. Unlike the fuzzy 
in the matrix model with a quintic term [60], the fuzzy CP^ in the present model 
is nonpcrtiubatively stable in the large-A^ limit despite the fact that it has larger free 
energy than the fuzzy in the same model. Thus the model provides a nonperturbative 
definition of a gauge theory on the fuzzy CP^. It would be interesting to investigate the 
field theoretical aspects of this model as in the case of noncommutative torus [41]. 

Prom the viewpoint of the dynamical generation of space-time, however, we should 
note that the fuzzy CP^ cannot be realized as the true vacuum since it has larger free 
energy than the fuzzy S^. This conclusion is in contrast to the results [53-55,57] obtained 
in the IIB matrix model [43], where four-dimensional space-time is shown to have smaller 
free energy than the space-time with other dimensionality. We should also note that the 
gauge group dynamically generated for the cases studied in this paper as well as in the 
previous work [59] turns out to be of rank one, although there is no reason for it a priori. 
Indeed ref. [78] presents an explicit model in which the gauge group with higher rank is 
realized in the true vacuum. In the IIB matrix model, we feel that supersymmetry plays 
an important role in obtaining four-dimensional space-time as well as a gauge group of 
sufficiently high rank that can accommodate the Standard Model. We would like to report 
on these issues in the near future. 
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A. Explicit form of the CP^ configuration 

In this section we present the explicit form of the representation matrix J"^™'"-* of the 
SU(3) algebra. This, in particular, provides us with the explicit form of the fuzzy CP^ 
configuration A^^ in eq. (2.7). 

For that purpose we need to introduce the so-called (anti-) symmetric tensor product. 
Let us denote the matrix element of the matrix A for the orthonormal states \i) and \j) as 
{A)ij = {i\A\j). The usual tensor product is defined by 

{ii,i2\A0B\ji,j2) = {ii\A\ji){i2\B\j2) , (A.l) 



- 11 - 



where |ji,i2) = |ii)b2)- The (anti-) symmetric tensor product axe defined through their 
matrix elements 

sym(^l ) ^2 |^®-B|ilj2)sym , (A.2) 
asym(«l,«2|^<8)-B|ilj2)asym , (A.3) 

where \ ji,j2)sym and \ji, j2)asym. are the orthonormal (anti-)symmetrized state defined by 

f |ji)|i2) for jl =j2 , 

[^i\jl)\j2) + \j2)\jl)) forii7^j2 , 



j2)sym 

\jl 1 J2)asym 



(A.4) 



V2 



(|ji)|j2) - |i2>|ji)) for jl / j2 , 



(A.5) 



respectively. The size of the matrices representing the symmetric tensor product {AfS>B)syni 
and the anti-symmetric tensor product (^(gii?)asym is k+iC2 and ^02 , respectively, where 
k is the size of the matrices A and B. The (anti-) symmetric tensor product can be gener- 
alized straightforwardly to a product of more than two matrices. Then the representation 
matrix of the (m, 0) representation is given as 



rp{m,0) 



sym 



(A.6) 



j-l m-j 

where denotes the fundamental (1,0) representation of the SU(3) Lie algebra, which is 
given explicitly as 



t4 




\ 


^1 


f^l o\ 


, h = 


0-10 


J 




[0 Oj 




^1 


\ 

-i 






[0 i 1 




The (0, n) representation can be obtained by simply replacing the fundamental represen- 
tation tfj, by the anti-fundamental representation, = — t* , in the above definition. 

In order to obtain the (m, n) representation, we have to define an orthonormal state 



) jm ) ^1 



) "'n/mix 



with a mixed symmetry such that it is symmetric with respect 



to the first m indices and the last n indices, separately, and anti-symmetric with respect 
to the exchange of one of the first m indices and one of the last n indices. (This symmetry 
is exactly the symmetry of the Young tableaux for the (m, n) representation.) We denote 



the tensor product defined with these states as (^i (g) • ■ • (g) (g) Bi 
the representation matrix of the (m, n) representation is given as 

m 

^ V I3 (g) • • • (g) I3 <^*M I3 <^ • • • <^ I3 ® I3 • • • 



-B„)mix- Then 



rp{rn,n) 



n 



i-1 



m-j 



k=l 



•^J mix 



(A.7) 



fc-i 



n—k 
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Note also that 



l3) + --- + (l3 



m terms 
+ (V ® V ® I3 «> • • 



13) + ••• + (13 



m(m— 1) terms 



n terms 



I3 <8) ® -^m) + ■ 



n(n— 1) terms 



l3) + ' 



2mn terms 

'4 1 1 

= ( + n) + -(m(m — 1) + n(n — 1)) + - x 2mn ) Ijv 



3' ' 3 

m{m + 3) + n(n + 3) + mn 



Ijv 



where we have used the formulae 



1 , 

3(13 



'-3jsym 



{tp <8) S^)asym = g (I3 <^ l3)asym , 
= {s,f = 1 13 . 



(A.8) 

(A.9) 
(A.IO) 
(A.ll) 



B. Perturbative expansion ciround fuzzy manifolds 

In this section we formulate the perturbation theory around the classical solution given 
either by (6.1) representing k coincident fuzzy CP^ or by (D.l) representing k coincident 
fuzzy S^. The single fuzzy CP^ and the single fuzzy are included as a special case k = 1. 

Let us evaluate the partition function Z = J dAe~^ around the classical solution 
Ap = Xp at the one-loop level. We define the measure of the path integral as 



8 N^-l 

dA=ii n 

/i=l a=l 



(B.l) 



where = J2a=i^ with being the generators of SU{N) normalized as tv{t"'t^) = 
Sab- We need to fix the gauge since there are fiat directions corresponding to the transfor- 
mation 



(B.2) 



where g is an element of the coset space H = \J(N)/\J(k). In order to remove the associated 
zero modes, we introduce the gauge fixing term and the corresponding ghost term 



Sg.f. = -ytr[X^,A^]2 , 
5ghost = -N tr:([Xf„ c][Ap,c]^ , 



(B.3) 

(B.4) 
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where c and c are the ghost and anti-ghost fields, respectively. These ghost fields take 
values in the tangent space of H. We perform the integration over ^4^ by decomposing it 
into the classical background and the fluctuation as = + The partition function 
can be rewritten as 

Z = YoliH) M J dAdc dc e-^t°t-i , (B.5) 
where the total action is defined by 

(B.6) 

and it is given explicitly as 

S'total = + S'kin + S'int , (B-7) 

^kin = ^iVtr +Nlv(c[Xx, [Xx,c]]) , (B.8) 

5int = -iVtr - ^ A^tr 

+1 i af^upN tr (i^^i.ip) + iVtr (c [X^, [A^, c]]) . (B.9) 

The normalization factor M = (27riV)-(^'-^')/2 in (B.5) can be obtained by following the 
usual gauge fixing procedure as in ref. [78]. The linear terms in cancel since is 
assumed to satisfy the classical equation of motion. Since the classical solution we 
are considering is proportional to a, we can rescale the matrices diS A^ a A^, c a c, 
c I— > ac, so that all the terms in the total action S'totai become proportional to a^. This 
means that the expansion parameter of the present perturbation theory is . The volume 
of the coset space H in (B.5), vol(if) = vol(U(A''))/vol(U(/c)), can be obtained by using 
the formula 

™'("«' = (p-l)!(;-2)....l!0! ■ 

We calculate the free energy W = — log Z as a perturbative expansion W = YljLo > 
where Wj = 0(q;^^^~-^^) comes from the j-loop contribution. The classical part is obtained 
as Wo = S[X], which is nothing but the action evaluated at the classical solution A^ = X^. 

Introducing the operator 

V^m"^' [X^,M] , (B.ll) 
which acts on a N x N traceless hermitian matrix M, the kinetic term (B.8) reads 

Skin = iVtr|i4(pA)'4 + c(^A)'c| . (B.12) 

The one-loop term can therefore be obtained as 

Wi = STr\og{N{Vxf} - log{vo\{H)Af} , (B.13) 
where the symbol Tr denotes the trace in the space of traceless hermitian matrices. 
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C. Perturbative calculations for the fuzzy CP^ 



In this section we focus on the case where the classical solution is taken to be the k 
coincident fuzzy CP^ (6.1). The results for the single fuzzy CP^ can be readily obtained 
by setting k = 1. 

C.l One- loop calculation of free energy 

Let us calculate the free energy up to one-loop. The classical part is obtained as 

Wo = -\N^a\n-l) , (C.l) 

N 1 

n=j = -{m + l){m + 2), (C.2) 



where we have defined 



and used the relation 

ffii/p fnvp' — 3 5pp' ■ (C.S) 

Next we evaluate the one-loop contribution Wi in (B.13). In order to solve the eigen- 
value problem of the operator {Px)^ , we introduce an analog of matrix spherical harmonics 
in the fuzzy case [59]. Let us denote it as {Ygt], where the indices s and t run over 
s = 0, 1, • • • , m and t = 1, • • • , (s + 1)^, respectively. It gives a complete basis for the space 
of n X n matrices. (Note, for instance, X^^o(* + -'^)^ ~ * given s, Ygt transforms as 

a (s,s)-type irreducible representation of SU(3) under the adjoint operation [T/^ , • ]. 
For more details of see refs. [28,62-65]. We also introduce /c x A; matrices e^"'*^ whose 
(a, h) element is 1 and all the other elements are zero. Then, as a complete basis ol N x N 
matrices, we define 

yif^ = Yst ® e(«'^) , (C.4) 

which satisfies the relation 

^^l^y{a,b)tyia\b')^^ = dss'Stt'Saa'Sw ■ (C.5) 

The eigenvalue problem of the operator (Vx)^ can be solved as 



{rx?ylf^=a's{s + 2)yifK (C.6) 



Note that J'oo'^'* for all the (a, b) blocks are the zero modes. In the k = 1 case, the zero 
mode should be excluded since are traceless. For k > 2 the tracelessness condition 
removes only one of the k"^ zero modes. Here we omit the rest of them by hand ^. Then 
the one-loop contribution Wi is given by 

m 

Wi = Zk^ ^(s + if log [Na^s{s + 2)] - log{vol(i?)A^} . (C.7) 



'^In fact these zero modes are responsible for the instabihty of the multi-fuzzy CP^ discussed in section 6. 
(See appendix D of ref. [59] for more in-depth discussions on this point in an analogous model.) However, 
the number of zero modes, fe^ — 1, is negligible compared with the dimension of the configuration space, 
which is of O(iV^). Indeed figure 6 shows that the results obtained by omitting the zero modes are in 
reasonable agreement with Monte Carlo results obtained before the decay of the multi-fuzzy CP^ actually 
takes place. 
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The one-loop free energy is obtained at large N as 

PFo + W^i~iV2 ^-^ + 61oga + log^^-^j , (C.8) 

where the rescaled parameter a is defined by (3.1). 
C.2 Derivation of the critical point 

In section 3 we observed a phase transition in Monte Carlo simulation starting from the 
fuzzy CP^ configuration. Wc can derive the critical point in the same way as in the 3d 
YMCS model [59]. Let us consider the effective action for a one-parameter family of 
configurations = /3r^"*'°^ 'Silk- At the one-loop level, it is obtained at large N as 

ri-ioop(^)^A^'|^(^-a^')+61og^ + log^-^| , (C.9) 

where P = PN*. The function of P on the right hand side has a local minimum for 

a > a^r^^'^ = ^k^ = 2.3094011 • • • x fci , (C.IO) 
v3 

which determines the (lower) critical point of the first-order phase transition. In fact it is 

known that the effective action is saturated at one loop in the large- limit in the case of 
fuzzy or fuzzy x [32,57]. This is the case also for the fuzzy CP^. Therefore, the 
critical point obtained above is expected to be correct to all orders in perturbation theory. 

C.3 One-loop calculation of various observables 

In this section we extend the perturbative calculation to various observables including those 
which are studied by Monte Carlo simulation in section 3 and 6. The zero modes, which 
appear for A; > 2, is omitted as in the evaluation of the free energy given in appendix C.l. 
We note that the number of loops in the relevant diagrams can be less than the order of 
^ in the perturbative expansion since we are expanding the theory around a nontrivial 
background. At the one-loop level, the only nontrivial task is to evaluate the tadpole 
ii^nhj) explicitly. 

C.3.1 Propagators and the tadpole 

The propagators for and the ghosts are given as 



ah s=\ t=l ^ ' \ / K( 

m (3+1)=^ . . 

(c)«(a)«)„ - E E E ^^^^ {y^% {yif") ^, • (c-^^) 

ab s=l t=l ^ ' \ / Kl 



where the symbol ( • )o refers to the expectation value calculated using the kinetic term 
Skin in (B.8) only. 
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Due to the symmetries, the tadpole (^^) can be expressed as 

with some coefficient c. Using the identity 

tr = ctv{X^X^) 

= -r-a m{m + 6) , 

the coefficient c can be determined by calculating the left hand side of (C.14). 
At the leading order in we have 



(C.13) 



(C.14) 



N 



- { triX^Af,) tr ( - i aUpaA^ApA^ 



(tr(X^i^)tr (c[X^,[i^,c 



(C.15) 



Using the fact that Xp_ is a linear combination of (y^=i j 1;;.), we can calculate (C.15) 
similarly to the previous section. After some algebra we arrive at 



tr (^X^(ij^)i-ioop) 



2^ 



(n^ - 1) . 



Using (C.14) we obtain 



loop — ' " - A„. ~ 



■x„ 



ri^mim + 3) ^ na^ ^ 



(C.16) 



(C.17) 



C.3.2 One-loop results for various observables 

Using the propagator and the tadpole obtained in the previous section, we can evaluate 
various observables easily at the one-loop level. 

The "extent of space-time" (-^tr(y4^)^) can be evaluated as 

1 



N 



1— loop 



1 

N 

a" 



ii{X^X^) + 2tr (x^(i;,)i-ioop) + (tr(i;.)2)o 

n^a^ ^ s(s + 2) 

s=l . 



1 . ox 4 - 1 
-m{m + 3) J ^ h 



At large with fixed a, we obtain 

1/1 



2a' 



(C.18) 



/iV\A^''^'^'^^Vi-ioop^ 3, ,2- 
The expectation value of the Chern-Simons term 

M = ^fp,ptT{A^A,Ap) 



(C.19) 



(C.20) 
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can be evaluated as 

6(n2 - 1 



—m(m + 3) + „ 
3 an^ 



t(M)i-.oop--^ + ^. 



(C.21) 
(C.22) 



At large A'' with fixed a, we get 

1 

m 

The observable (itr(F^,,)2) can be calculated in a similar manner, but it is easier to make 
use of the Schwinger-Dyson equation 



ltr(F,.)2 + 3aM^ = ^ (l " 



from which we obtain 
1 

N , i_ioop 

Combining (C.21) and (C.24), we get 
1 



= 8(1-^1 -3a(M)i_ioop 



2a^ 



- 10 



1— loop 



"6fc ^ 2 



(C.23) 



(C.24) 



(C.25) 



C.3.3 An alternative derivation 



Since tr(i<'^^)^ and M are the quantities that appear in the action 5, we can obtain their 
expectation values easily by using the free energy (C.8) calculated for the k coincident 
fuzzy CP^. Let us consider the action 



S{l3uh;a) = Ntv (^-^[Ai„A,f^ + (32N^aM , 



(C.26) 



where we have introduced two free parameters /?i and (32, and define the corresponding 
free energy by 



By rescaling the integration variables as A^, we find 

W{Pi,(32; a) = 2{N^ - 1) log/?i + W{1, 1; a(32Pih ■ 



(C.27) 



(C.28) 



Then (tr(F/,,,)2), (M) and (S) can be obtained by 



(M) 



d(3i 

1 dWiPi,P2;a 



3a dW{l,l;a) 



1 / 1 



802 



1 dW{l,l-a) 



da 



da 



1 



a dW{l,l;a) 



A^2 / 4iv2 



da 



(C.29) 
(C.30) 

(C.31) 
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Using the one-loop result 

W{1, 1; a)i_ioop = --N^a\n - 1) + 3k^ ^(s + l)^ log [Nah{s + 2)] - log{Yo\{H)Af} , 

(C.32) 



6 

s=l 



which follows from (C.l) and (C.7), we can reproduce (C.21) and (C.24). 
C.4 All order results from one-loop calculation 

Taking advantage of the one-loop saturation of the effective action mentioned at the end of 
appendix C.2, we can obtain all order results for various quantities in the large- limit by 
simply shifting the center of expansion in the one-loop calculation [37] . Similar calculations 
have been done previously in the 3d YMCS model [66]. 

Since the free energy and the effective action are related to each other by the Legendre 
transformation, we can obtain the free energy by evaluating the effective action at its 
extremum. We consider the expansion around a configuration = [3T^ Ifc- The 
value of /9 that gives the local minimum of the effective action can be obtained by solving 

^ri-ioop _ ^2 [2 _ -^2) + 11 = 0. (C.33) 



d(5 yV"^ ^ ' /3 

The solution exists for a > a^^^ ^ = -^k*, and it is given explicitly as 



P = f{a)^jU + Vl + 6 + J2-6 + -^==j, (C.34) 



where 

S = a"t(96fc)5 < 



. 256A;\" I I 256k \ . ,^ 



9a^ j \ V 9a^ 
At large a, the solution (C.34) can be expanded as 

/(a) = a I l-^c,a n = « - _ - — J. (C.36) 

Plugging this solution into (C.9), we obtain the free energy to all orders as 

It.. , 9 9k 63P 1485P 

_„,.__ + 61oga + log— -j-^ — (C.37) 

Using (C.29), (C.30) and (C.31), we obtain the all order results 

/ 1 , ,n\ 2a ^, ,o 2oL^ 108A; \h\2k^ 26730A;3 

(C.38) 

1 Viptf 2c? 6 36A: 504A;2 SgiOA:^ ^ 

/i^i 3k 3k a ou^ a^'^ 

^ ln^ 1 n3 7 9A; 126^^ 4455^^ 

^(S> = +2 = -_ + - + _ + _ + _ + ... .(C.40) 
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Similarly to the case of the 3d YMCS model [66], we can also calculate various ob- 
servables directly to all orders in perturbation theory in the large-N limit. For instance, 
the observable -^^■i^tr(A^)^^ is obtained at one loop as (C.19), whose first and second 
terms correspond to the classical and one-loop contributions, respectively. As we see from 
(C.18), however, the one-loop contribution comes from the tadpole diagram, which is not 
one-particle irreducible. Therefore, by replacing a by f{ot) in the classical contribution, 
we obtain the all order result as 

1 /I , , ,o\ 2f(af 2c? 4 30/c 432A;2 7722A;3 

/ x,{A^f\^J^ = — - — . (C.41) 



Viv\A^ / 3A; 3fe 6? a^o a^^ 

D. Perturbative calculations for the fuzzy 

In this section we perform the perturbative analysis for the k coincident solution 



X = ^(fcs2) ^ f a^A*"^ <^ Ifc for /x = 1,2,3, 
I otherwise , 



which generalizes (2.11). The total size of the matrix is now given by AT = nk. Let 
us calculate the free energy as a perturbative expansion W = Yl'jLo^j^ where Wj = 
0(a^(^'~-')) comes from the j-loop contribution. Prom (B.13) the free energy is obtained at 
the one-loop level as 

Wo + Wi = —{n^ - 1) + 3k^ + 1) \og[NaH{l + 1)] - log{vol(ff)AA} 

1=1 

"^'(-2|^ + ^^°^^ + ^°S^-t) ' 

where the rescaled parameter a is defined by (4.1). 

Similarly, the effective action for a one-parameter family of configurations 

^ ^ f/34"^®lfc for ^ = 1,2,3, ^^ ^^ 

^ 1 otherwise 

can be obtained at the one-loop level as 

r,-W/5) = iv'^(iJ-|^ + 6iog,5 + iogf -fj , (D.4) 

where (3 = (3N^ . The effective action has a local minimum 



/3 = g{a) = ^ [ 1 + + J 2 - £ + -^=== 



~ 4, ,,2^1 f/^ I 1024)^2 \^ / I 1024fc2 \ 3 \ 
e = a-.(384/c^). Hl + ^l-^j + ^1 - ^1 - , (D.5) 
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if a is larger than the critical point 



a(f = = 3.2659863 ■ ■ ■ x . (D.6) 

Since the effective action is saturated at the one- loop level in the large- limit [32], the 
critical point obtained above should be correct to all orders in perturbation theory. 

The propagators for and the ghosts are exactly the same as (C.l) and (C.2) in 
ref . [59] . The tadpole is given by 

(^m)i-Ioop = ~Jp:^^i^ ■ ^^■'^) 
The one-loop results for various observables are obtained as 

— M i_,oop = ^ ^ + 1 - — 

^tr(F^,)2^^ ^ = 8 (^1 - ^) - 3a(M)i_ioop 



= \a\n' - 1) - 10 + ^{8{3k' - 1) - 6k'} 

^(5)i-ioop = I (^^tv{F^u?^^ ^ + a(M)i_ioop 

7 



24/c2 2 



1/1 ,2\ 2|1,2 6fc2(n2-l) , 8 "^21 + 1 



N \N ' /i_ioop l4A^ ^^"^ iVra2a4 ^ /(i + 1) 



-2/1 6 16 , 

Exploiting the one- loop saturation of the effective action in the large- limit, we can 
calculate various quantities to all orders. The free energy is obtained as 

1 ~ 1 15 36^2 lOOSA;^ A7520k^ 
iV^^"-24F + '^°S« + log^-^-^^ ^3 ^T2 

(D.12) 
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by plugging (D.5) into (D.4). Various observables are calculated as 

1 1 . ,s ^ oi^ 7 36/c2 2016fc^ 142560/0*^ 

(D.13) 

1/1 . . .2\ 1 /~^2_ 6 180A;2 10368A;4 741312fc6 



tr(^M)7 - ^^(«)^ = 4^ - i - - - - • • • ' (^-1^) 

1 1 . 6 144fc2 8064A;'^ 570240^*^ , 



(^tr(F,.)^ 



6fe2 6fc2 a ct' cfi o^'^ 
1 _ „ 432A;2 24192^"^ 1710720^6 
+8=2p-^'-^i 5^ ^ 



(D.16) 
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